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Is Mira a magneto-dusty rotator? 
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ABSTRACT 

We investigate the possibility that a magnetic field may be present in the star o— Ceti 
(hereafter, Mira) and that the field plays a role in the star's mass loss. The model 
presented here is an application of an earlier derived theory that has been success- 
fully employed for intermediate and high-mass evolved stars, and is now extended 
to the low-mass end. The modelling shows that it is possible to obtain a hybrid 
magnctohydrodynamic-dust-driven wind scenario for Mira, in which the role of a mag- 
netic field in the equatorial plane of the star is dynamically important for producing a 
stellar wind. The wind velocity and the temperatures obtained from the model appear 
consistent with findings elsewhere. 

Key words: MHD-stars: AGB and post-AGB - stars: winds, outflows. 



1 INTRODUCTION 

Mira (omicron-Ceti) is a relatively close evolved star that 
is in its asymptotic giant branch (AGB) infancy. It is the 
template upon which the class of stars called Mira vari- 
ables are based, which are characterised most strikingly, by 
large-amplitude long period v ariability. Mira va ries with a 
period of about 332 days (e.g. lLopez et al]|l997f ). Observa- 
tions of emission due to vibrational-rotational transitions in 
the CO molecule have not only revealed that the envelope 
of Mira is mildly asymmetric, possibly due to a bipolar out- 
flo w, perhaps due to an eq uatorial magnetic field (see for 
e.g. iPlanesas et ai1ll990al lbl; lLopez et alii 19971 ) but also CO 
transition line observations have also enabled a variety of dif- 
fere nt estimates of the wind velocity ranging from 4.8 km/s 
( see l Young! 19951 ) to 2.5 km/s for the envel ope's expansion at 
intermediate distance s of 100— 1000-Rq (see lRvde et aljbooot 
iRvde &: Sc hoier 2001 ). Elsewhere , there is evidence for dou- 
ble winds fsee lKnapp et al.lll998l ) with a fast wind of about 
6.7 km/s, and a slow component of about 2.5 km/s. Mean- 
while, observations in the infrared at around 11/im, have 
revealed that the inner rad ius for the dust around Mira is 
at ar ound 3-Ro (see for e.g. iBester et al.lll99ll ; lLopez et al.l 
1 19971 ). where R® is the photospheric radius. T hese observa- 
tions have also been modelled with success (see lLopez et al.l 
1 19971 ) using an additional dust shell at 12J?o, as well as with 
dust clumps with photospheric spots. 

SiO-maser observations of the molecular shell close to 
the photosphere of Mira have revealed that the star may 
be rotating with a period of about 89 x sin(i) years (see 
ICotton et al.l I2006T ) . where i is the angle between the line 
of sight and the rotation axis. SiO maser polarisation stud- 
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ies on the other hand, indicate that Mira may harbour a 
predominantly radial ma gnetic field in its atmosphere (e.g. 
ICotton et al]|2004l2006l ). Elsewhere, recent observations of 
masers have been instrumental in establishing the dynamic 
role that magneti c fields play in the atmospheres of evolved 
stars (see for e.g. [ Amiri et al.ll2012l ; IVlemmings et all 1201 ll; 
Herpin et al.l 12009 ) and theoretical efforts (see Busso et al.l 
20071) have delineated their role in transport of material in 
AGB stellar interiors with surface fields of ;S 20 G. 



Imaging of Mira in the near infrared (IR) and optical 
wavelengths has provi ded tools for measuring the diameter 
of the star (see for e.g. Haniff et al.lll995l ; lvan Leeuwen et al.l 
ll997l) . |Perrin et alj (|2004l ) find that by varying the opacity of 
the molecular layer just ahead of the photosphere, they were 
able to account for the apparent changes in the diameter of 
Mira, arriving at an estimate of w 35Oi?0, over the entire 
variability cycle of Mira. 

From the brief discussion of the literature presented 
above, it appears that the role that magnetic fields may 
play in shaping the outflow and influencing certain dynamic 
features of the envelope cannot be ignored. In the current 
study we present a rudimentary model integrating the ef- 
fects of rotation, an equatorial magnetic field and the usual 
dust-driving picture into one cohesive scenario for Mira's 
outflow. The work prese nted here represents an extens ion 
of our earlier theory (see iThirumalai fc Hev3l20ld . \20lj ) to 
the low-mas s end of AGB stars ; viz., Mira which is about 
1.5M Q (e.g. iMartin et al.l I2007T ). Now that stellar rotation 
may have been detected in Mira, the aim of the current 
work is to raise the question, how important is the magnetic 
field in the stellar outflow and can Mira be a magneto-dusty 
rotator, given the current observations? 
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2 THE HYBRID- WIND MODEL 

We confine our attention to an axisymmetric model in the 
equatorial plane of the star. Therein, the velocity fields of 
the gas and the dust as well as the magnetic field are as- 
sumed to have radial and azimuthal components. However, 
the azimuthal components are functions of purely the radial 
distance from the star; this is the c entral assumption be- 
hind the canonical Weber-Davis (see lWeber fc Davislll967l . 
hereafter WD) model for our Sun. There are two fluids; the 
gas which carries the magnetic field and secondly the dust, 
which moves through the gas, dragging the gas with it. The 
two fluids are coupled by a drag term in the usual way. 
The dust-to-gas ra tio is kept small around ~ 1/355 (e.g. 
iDanchi et al.l Tl994T ) and the individual dust grains are as- 
sumed to be spherical in shape. The gas is assumed to have 
a polytropi c equation of state, with a poly tropic exponent 
7 > 1 (see IThirumalai fc Hevlll2010l . 12012 ). where a value 
of unity represents the isothermal limit. We assume ideal 
MHD where there is no Lorentz force in the fluid. The input 
parameters for model of Mira's wind are listed in Table 1. 
A steady-state description of the gas velocity in the hybrid 
wind in the equatorial pl ane of Mira can be written as (see 
IThirumalai fc Hevllboiol ). 



(1) 



dw w N(w, x) 
dx x D(w, x) ' 

where, w — u/ua is the gas speed normalised using the 
Alfven speed and x — r/rA, is the radial distance expressed 
in units of the Alfven radius. Hereafter, the subscript 'A' 
refers to values of the different variables at the Alfven radius. 
The quantities N(w, x) and D(w, x) are the numerator and 
denominator respectively and are given by, 
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along with 7 uniquely determine the 

locations of the critical points, and hence the morphology of 
the family of solutions of Eq. fTJ. Here Ta is the gas tem- 
perature at the Alfven radius, k is the Boltzmann constant 
while m v is the mass of a proton. The presence of the Heav- 
iside function in Eq. represents the formation of dust at 
the loc ation x = xg- The dust ve locity profile is then given 
by (see IThirumalai fc Hevlll2010l ). 



v(r) — u(r) + 



+ 4 



— 1 1 



2 
O-th 



1/2 



V 



(4) 



/ 



where a t u is the thermal speed given by a t h = \/2kT/ fim u 
and firriu is the mean molecular mass of the gas and rid is 
the dust grain number density, which is assumed to be given 
by, n d m,d/ p ~ (8), with (5), the average dust-to-gas ratio in 



Table 1. Various parameters for modelling Mira. The variable 
parameters are listed as such. 



Parameter 


Symbol 


Value / Comment 


Mass 


M* 


~ 1.5M 


Radius 


Ro 


- 464R Q 


Mass loss rate 


M 


~ 3 X lCT 7 M /yr 


Surface magnetic 






field strength 


B 


variable 


Bulk surface gas 






velocity (radial) 


"0 


variable 


Surface temperature (effective) 


To 


~ 2500K 


Stellar rotation rate 


n 


2tt/ [89 X sin(j) yr] 


Rotation axis angle 


i 


variable 


Surface escape velocity 


Vesc,0 


3.52 X 10 6 cm/s 


Polytropic exponent 


7 


> 1 


Alfven radius 


TA 


variable 


Alfven speed 


"A 


variable 


Dust Parameter 


r d 


variable 


Dust grain radius 


a 


spherical grains 


Dust grain mass 




~ A/3na 3 p d 



the wind. Eq. fl]) is solved num erically on a computer (see 
IThirumalai fc Hevlll2oT(il . l2012h . In the following section the 
results are presented alongside a discussion. 



3 RESULTS AND DISCUSSION 

Along with the basic ingredients of the model listed in Ta- 
ble 1, the set of parameters {-Bo, uo, u a , r a , 7, Fd, i} are var- 
ied until a critical solution to Eq. ([T]) is obtained that satis- 
fies the following criteria. 

1. The solution passes through all three critical points; 
the sonic point, the radial Alfven point and the fast point 

2. The solution is continuous through the radial Alfven 
point. 

3. The velocity profile starts at the base of the wind sub- 
sonic and attains a super-Alfvenic terminal velocity at large 
distances 

4. The temperature range in the dust condensation region 
(within a few stellar radii from the photosphere) is consistent 
with observations. 

5. The gas terminal velocity is consistent with observa- 
tions. 

This optimisation procedure carried out in tandem with 
integrating the differential equation in Eq. (O results in a 
picture of a hybrid MHD-dust-driven wind for Mira as shown 
in Figure [1] The critical wind solution is shown as the solid 
red line, which comprises of two parts, L\ and C2. These 
two lines intersect at the dust formation radius r d , which in 
the current model is located at 3Ro- In the region r ^ r d 
Eq. |T} is integrated with the presence of the dust parame- 
ter T d , which represents dust grain drag acting on the gas; 
it is this part of the solution that is labelled C2. On the 
other hand, inside r ^ r d there is no dust and as a result 
Eq. |T} is integrated in this latter region without the dust 
parameter, as a pure WD solution. This part of the solu- 
tion is labelled as L\. Thus in the region Ro ^ r ^ r^, 
a pure WD mechanism is responsible for transport of stel- 
lar material. Together, L\ + C2 forms the combined hybrid 
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Figure 1. Hybrid wind solution is shown for Mira with param- 
eters ua ~ 0.14i) eS c,0i T A ~ 43.47-Ro and for S3 0.06 and 
remaining parameters as given in Tables 1 and 2. The red solid 
line (Li + C2) traces the hybrid MHD-dust-driven wind solution 
for Mira. The decreasing blue solid line traces the temperature 
and should be interpreted using the right hand y-axis. 



MHD-dust-driven wind. The radial Alfven point and the 
fast point lie nearly coincident upon one another; shown by 
the black dashed vertical line at r « 43i?o- Theoretically 
speaking, if we assume that dust condensation occurs at the 
photosphere, then we would obtain the combined solution 
Ci + C2 by integrating Eq. |T} with the dust parameter, 
over the entire domain Rq r ^ 100i?o • However, in the at- 
mosphere of Mira, dust formation occurs at the dust conden- 
sation radius r — r d and not at the photosphere. Therefore, 
inside the dust radius, r ^ r d , the wind starts off at the pho- 
tosphere at some velocity uo, and proceeds outwards along 
the trajectory L\, and after dust condensation at r = rd, the 
solution then switches to proceed outward along C%, rather 
than continue along the unphysical solution L2. 

Thus, it is seen explicitly, that without the onset of 
dust formation at r — rd, there would not be any efflux 
from Mira, since there is only one solution (C2) that passes 
through all three critical points, and emerges super- Alfvenic 
at large distances. As can be seen, the gas terminal velocity 
is about 5 km/s in the equatorial plane, which is in reason- 
ably good agree ment with estimates of the wind velocity of 
Mira fe.g. lYoundlT995l ; lMartin et alj|2007h . 

The dust velocity profile is also shown in Figure [T] as 
the brown dot-dashed line. This is computed according to 
Eq. @, once the gas velocity profile is known. The dust ve- 
locity is slightly larger than the gas velocity, as expected, 
in a dust-driven wind. The red dashed line represent the 
azimuthal velocity profile of the gas. This profile is typical 
for a magneto-centrifugal wind. The blue solid line shows the 
temperature profile in the atmosphere of Mira and should be 



interpreted using the right hand side axis. The photospheric 
temperature is about 2500 K for the model shown. Finally, 
the green shaded area shows the so called hybrid region. This 
region is bounded to the left by r = r d n ~ 2.6 7?o, this is the 
lower limit for the inner dust radius as given bv lDanchi et al.l 
(|1994| ). To the right, the hybrid region is bounded by the 
sonic point of the pure WD magneto-centrifugal model. The 
sonic point of the hybrid model lies at r ~ 5.61i?o; just 
inside the hybrid region. As seen in our earlier study (see 
iThirumalai fc Hevll I2OT0I ) , one of the ways in which a hy- 
brid wind is possible is if the sonic point of the hybrid 
model lies within the sonic point of the pure WD model; 
i.e., r s < rY D ■ Moreover, the dust formation radius must 
also then lie inside r ^ r s . The temperature in the hybrid 
region can be inferred from the temperature profile shown, 
to be about 1000 5 T < 1300 K, which is well within the 
observed range of temperatures at this distance from the 
star ( e.g. lLopez et aTlll997l : iBester et al]|l99ll ; iDanchi et~ai1 
ll994l ; |Perrin et al.ll20o3 v 

Finally the optimised values for the different variable 
parameters in the model are given in Table 2. Notice that the 
surface radial magnetic field at the photosphere is obtained 
to be about 4 G, which is w i thin t he ra nge of field st r engths 
estimated by iHerpin et alJ (|2006h and iBusso et all (|2007f ) 
for AGB stars. 

Presently, we turn our attention to the question of hot 
spots on the photosphere of Mira, and the related question 
of the influence of the spot on the stellar wind, ahead of the 
spot in the atmosphere. Figure [2ja) shows a hybid MHD- 
dust-driven wind model where the photospheric tempera- 
ture is increased to about 2700 K. However, in this second 
model instead of formulating a hybrid wind model as be- 
fore, with dust formation occurring within the sonic point, 
we chose a different scenario. Our motivation here was to in- 
vestigate the possibility regarding dust form ation at around 
127?o, at which distance lLopez et alJ (|l997T ) model a second 
dust shell. Within the framework of the current theory, we 
found that the only way to have a viable wind scenario with 
dust formation at this distance, was by collocating the dust 
formation radius with the radial Alfven point. However, in 
this case, as can be seen from Figure [2] the solution is not 
smooth at the radial Alfven point. In this scenario, the solu- 
tion need not be smooth at r = ta, as this also happens to 
be the dust formation radius, hence a kink in the solution 
is allowed. It is possible that the dust formation can occur 
outside the fast point, by fine tuning the parameters of the 
model. This would place the dust formation radius in the 
vicinity of about 20 — 307?o- While we do not show such a 
scenario in Figure [5] it is to be acknowledged that it may 
be possible within the framework of the hybrid MHD-dust- 
driven wind theory. Our objective here was to formulate a 
model with dust formation occurring at a distance of about 
127?o, to be co nsistent with othe r models for dusty shells 
or clumps (e.g. lLopez et al] [l997). while obtaining temper- 
atures in agreement with observations. 

In Figure [T] it can be seen that the temperature at the 
left bound of the hybrid region (r = r d n = 2.6i?o) is about 
~ 1300 K. This same temperature is calculated to be at 
a distance of about 2.9Ro in Figure [2] Thus, ahead of a 
hot spot in the atmosphere, the viable dust formation re- 
gion moves further out from the photosphere, as would be 
expected. In Figure [2] the temperature at the sonic point 



© 2011 RAS, MNRAS 000, \T§5\ 



4 Anand Thirumalai and Jeremy S. Heyl 




r/Ro r/Ro 



Figure 2. (a) Hybrid wind solution is shown for a scenario with a 
hot spot at the photosphere of Mira. This model has parameters 
ua ~ 0.15i) esc ,0i i"A ~ 29.21i?o and ss 0.1 and remaining 
parameters as given in Table 1. The surface magnetic field was 
found to be Bo Ri 1.15G. The decreasing blue solid line traces 
the temperature and should be interpreted using the right hand 
y-axis. (b) Gas velocities obtained by perturbing the solution in 
the vicinity of r \ and (c) the corresponding dust velocities. 



is seen to be about ~ 900 K. Thus the range of tempera- 
ture in this model is similar to that for the hybrid region 
in Figure [1] therefore it is still possible to formulate a hy- 
brid wind scenario with dust formation occurring within the 
sonic point, as was shown in Figure [T] However, as men- 
tioned earlier, in the current model our goal is to locate a 
seco nd dust shell out side the sonic point at around 12i?o 
fsee lLopez et al.lll997l ) and not within it. 

In this scenario, with a short-lived hot spot on the pho- 
tosphere, the temperature in the gas is about 800 K at the 
radial Alfven point (r = ta)- Thus, the observed tempera- 
ture for t he dust of about 60 K at this distance is quite fea- 
sible fe.g. lLopez et al.1ll997l ). Overall, for the lifetime of the 
hot spot, it is possible to sustain a hybrid wind with dust 
forming at 12iio. A pure WD mechanism transports stel- 
lar material from the photosphere, through the sonic point 
and out to the radial Alfven point. After dust formation at 
this location, the MHD-dust-driven wind then negotiates its 
way through the remaining critical point; the fast point and 
leaves the star super- Alfvenic. Comparing Figure [2{ a) with 
Figure [1] reveals that the hybrid model is quite sensitive to 
changes in the photospheric temperature and the magnetic 
field. For the model shown in Figure [5] the surface radial 
magnetic field strength was found to be w 1.15 G. The dust 
parameter (r^) was found to be about 0.1; a little higher 
than in the previous model (c.f. Table 2). With these pa- 
rameters, the terminal wind velocity is also concomitantly 
higher, at around 11 km/s in the equatorial plane. In corn- 



Table 2. Optimised values of the variable parameters for the star 
Mira. 



Parameter Value 



B S3 4.11 G 

u « 5.89 X 10" 6 km/s 

ua 4.93 km/s 

r A & 43.47R 

7 « 1.06 

Y d w 0.06 

i 26.60° K, 0.46 rad 



parison to the time for one stellar rotation, since the hot-spot 
at the photosphere would be short-lived, therefore the hot 
spot model described above would be valid for the duration 
of the spot. 

In this regard, there are some concerns that are inher- 
ited from collocating the dust formation and radial Alfven 
points. In the main, dust formation and growth is a stochas- 
tic process that is not completely understood. Particularly, 
for the wind ahead of a short-lived hot spot in the atmo- 
sphere, the physical processes may be quite dynamic over 
time-scales comparable to the spot lifetime. As a result small 
changes in the wind velocity in the vicinity of the radial 
Alfven point can result in a drastically different wind veloc- 
ity profile and indeed the wind may not be able to navigate 
through the fast point. This is shown in Figure E^b) and 
(c). We have perturbed the gas velocity by a small amount 
0.01%) ahead of the radial Alfven point and then inte- 
grated the perturbed solution outwards. Figure [^b) shows 
the unperturbed solution as the red solid line and the per- 
turbed wind solutions as the dashed lines. We see that this 
perturbation of the gas velocity results in the wind becom- 
ing either a failed wind, where it does not pass through the 
fast point (the dashed lines below the solid red line) or the 
wind becomes an unphysical double valued wind (dashed 
lines above the the red solid line that form loops). In either 
case, this suggests that minor perturbations of the wind ve- 
locity around the radial Alfven point will drastically change 
the nature of the outflow. Figure [2jc) shows the correspond- 
ing dust velocity profile for both the unperturbed and the 
perturbed solutions. Should the dust follow any of the per- 
turbed trajectories then it will not leave the star and can 
lead to instabilities in the flow. Therefore, while it is possi- 
ble that in the atmosphere of Mira, conditions ahead of a 
mild hot spot may result in a hybrid wind with dust for- 
mation at around 12i?o, these conditions would be limited 
by the life-time of the spot, and as such cannot be clas- 
sified as steady-state, even when compared to the stellar 
rotation time. Such a localised phenomenon ahead of a hot 
spot may be able to account for dumpiness of the dust dis- 
tribu tion and even spa tially distinct dust shells around Mira 
(e.g. lLopez et al.lll997l '). We shall end this discussion by cau- 
tioning the reader that fully dynamic 2- or 3-D modelling 
would be needed for establishing the importance of such dy- 
namic short time-scale phenomena, which cannot be entirely 
gleaned from the simplistic theory presented here. The con- 
clusions of the current work are summarised in the following 
section. 
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4 CONCLUSION 

In this letter we have modelled the stellar wind of Mira as 
a hybrid MHD-dust-driven wind. The study was motivated 
largely by hints of the discovery of a magnetic field in the 
star in addition to a possible detection of stellar rotation. 

Since the surface magnetic field strength as well as the 
rotation rate are not exactly known, these along with a few 
other physical parameters such as the dust parameter (Td) 
were treated as variable parameters that were fine-tuned, in 
order to arrive at a suitable stellar wind model for Mira. We 
obtained a stellar wind with a terminal velocity of about 

5 km/s in the equatorial plane. This value for the wind ve- 
locity is in agreement with current estimates. The stellar 
effective temperature was taken to be about 2500 K and we 
obtained a surface magnetic field of about ~ 4 G, in good 
agreement with current estimates. We also arrived at an es- 
timate for a hybrid region in the inner atmosphere of Mira. 
This region depicts the range of distances at which dust 
condensation can occur in order to produce a hybrid wind. 
The temperature profile obtained for this hybrid region is 
consistent with models of the atmosphere of Mira. 

For the purpose of modelling dust shells around Mira 
at several stellar radii, we employed a model of a hot-spot 
that alters the wind dynamics ahead of the spot in the at- 
mosphere. With this we were able to show that it may be 
possible to facilitate dust formation at around 12i?o by col- 
locating the dust formation and radial Alfven points, result- 
ing in a hybrid wind. However, because such a spot may be 
short-lived the dust and gas velocities in the vicinity of the 
radial Alfven point are likely to change over short timescales, 
resulting in failed wind solutions and possibly dynamic in- 
stabilities. This may form form localised regions of dusty 
clumps, or even distinct shells with slightly enhanced dust 
density. 

It is however to be remembered, that the current model 
is an idealised steady-state description of a more com- 
plicated and dynamic problem. As such, the results ob- 
tained here only convey the possibility, that Mira may be 
a magneto-dusty rotator. Dynamic multi-D modelling of its 
wind is required for obtaining ultimately, a truer description 
of the nature of its efflux and of its atmosphere. 
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